SCALING LIMITS AND EXIT LAW FOR MULTISCALE 

DIFFUSIONS 
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Abstract. In this paper we study the fluctuations from the limiting behavior 
of small noise random perturbations of diffusions with multiple scales. The 
result is then applied to the exit problem for multiscalc diffusions, deriving 
the limiting law of the joint distribution of the exit time and exit location. 
We apply our results to the first order Langevin equation in a rough potential, 
studying both fluctuations around the typical behavior and the conditional 
limiting exit law, conditional on the rare event of going against the underlying 
deterministic flow. 



1. Introduction 

Let T > be given and consider a small random perturbation of dynamical 
system by a Wiener process. In particular, consider the d-dimensional process 
X e = {Xt, < t < T} satisfying the stochastic differential equation (SDE) 

(1.1) dX e t =b e (Xt)dt + ^a e (Xt)dW t , X^ = x , 

where e J, and Wt is a standard d-dimensional Wiener process. The functions 
b e {x) , a € {x) are assumed to be sufficiently smooth. 

If b e (x) — > b{x) and a e (x) — > cr(x) as e | 0, where b(x) and <j{x) are nice functions, 
then asymptotic behavior such as law of large numbers, central limit theorems and 
large deviations have been extensively studied in the literature, e.g., [H[T2] and the 
references therein. Scaling limits of (jl.ip under the effect of different perturbations 
of the dynamics and of the initial condition are also studied in the recent article 

&■ 

In this article, we assume that the functions b e (x) and <J € (x) are fast oscillating, 
in particular we set b e (x) = |6 (x, f ) + c (x, |) and cr e (x) = a (x, |), where 5 = 
<5(e) I as e J, 0. The functions b(x, y), c(x, y) and a(x, y) are assumed to be smooth 
and periodic with period p in every direction with respect to the second variable. 
Homogenization of such equations has been studied extensively in the literature, see 
for example [3J [53] . Large deviations were studied in [TU] and related importance 
sampling schemes were developed in [5J[7]- Moreover, special cases of this general 
equation (e.g., with b(x,y) = —WQ(y),c(x,y) = —W(x) and a(x,y) — constant) 
have been suggested as models for studying rough energy landscapes that describe 
certain proteins and their folding and binding properties. A representative, but by 
no means complete, list of references is [H [71 [HJ [TH1 [Ml US- 

Our goal in this paper is twofold. First, we study scaling limits under different 
perturbations of the drift and of the initial condition. We are interested in fluctua- 
tions around the typical behavior of X\ as e, S J, when both the initial condition 
and the drift follow a scaling limit in finite time. It turns out that depending on 
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the scaling and on the order that e and 8 go to zero, we have different limiting 
behavior. The result is presented in Theorem 13. II It is interesting to note that, in 
contrast to the case without fast oscillations, in the case considered here, additional 
drift terms may appear in the equation that the fluctuation process satisfies, see 
Remark 13.31 At this point we mention the articles |9j [TUl Q31 [17] for some related 
moderate and large deviations results, even though the fluctuations analysis done 
in the current paper is not covered, as far as the authors know, by the existing 
literature. Theorem 13.11 allows us then to study the exit distribution in the limit as 
e, 6 I (see Theorem l4.1|) in the case in which the typical behavior of X e exits the 
domain transversally in finite time. 

Second, we are interested on the effect that perturbations by small but fast 
oscillations of small noise dynamical system have on exit time for such diffusions 
conditioned on rare events, see Theorem 15.31 We investigate this question in the 
case of the first order Langevin equation for both a periodic and for a random 
rough potential, see Remark l5.4l It turns out that the limiting law of the exit time 
conditioned on the event of going against the deterministic flow, appropriately 
normalized, follows Gaussian distribution with enhanced variance (as compared to 
the small noise not oscillating case) due to the fast oscillations, sec Remark 15.51 

The rest of the paper is organized as follows. In Section [2] we establish notation 
and mention examples and preliminary results that will be used throughout. Section 
|3] contains the corresponding central limit theorem, whereas Section |4] contains the 
analysis of the joint limiting law for the exit time and exit point. In Section [5] 
we apply the results of Sections [3] and [5] to the first order Langevin equation in a 
rough environment. In particular, we state the related central limit theorem and 
study the conditional exit law of a one dimensional small noise diffusion process in 
a rough environment in the limit as the fluctuations and noise intensity go to zero. 



2. The set-up 



Let T > be given and consider the d-dimcnsional process X e 
satisfying the stochastic differential equation (SDE) 
(2.1) 



{X t e ,0 <t<T} 



dX! 



b Xf 



XI 



xt 



Xi 



Xi 



dt+sfeo I X 



X', 



dWt. 



with initial condition given by X* = x + e a2 ' 2 ^. Here C is a family of random 
variables that converges in distribution to £° as e — > 0, 5 = 5(e) — > as e — > and 
Wt is a standard d-dimensional Wiener process. Also, we assume that the functions 
b, c, *i> e and a satisfy the following conditions: 

Condition 2.1. (i) The functions b(x,y),c(x,y),o~(x,y), and ^^(x^y) are, 

for each e > 0, periodic with period p in the second variable, C (y) in y and 
C 2 (EL d ) in x with all partial derivatives continuous and globally bounded in 
both variables. Here y — T d denotes the d- dimensional torus, 
(ii) As e — > v I /£ — > ^ uniformly in each variable and satisfies the same 

regularity conditions as any ty e . 
(Hi) The diffusion matrix o~a T is uniformly nondegenerate. 
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We are interested in the following cases of interaction 

(2.2) lina| = (°° ***** 

eio (5 (^7 g (0,oo) Regime 2, 

Here 7 is taken to be 7 = 00 in Regime 1. 

We borrow some notation from [5], where the large deviations principle for 
SDE (|2.1[) was established, in order to present our results. 

Definition 2.2. For each one of the Regimes i = 1, 2 defined in i2.2\) . and x G M d ; 

define the operators 

Cl = b(x, •) ■ V v + -fr [<r(x, >(a;, -) T Vy] , and 

4 = [ 7 6(x,-) + c(x,-)] ■V a +7^r[ ( x( a; ,-)a(x,-) T V2] . 

-For eac/i x G R d , t/ie domain of C x is given by T>(C X ) = C 2 (y), for i = 1,2. 

Note that the existence of a unique smooth invariant measure for the operator 
C x , i = 1,2, is immediately implied by Condition 12.11 Wc impose the following 
condition for the invariant measure in Regime 1: 

Condition 2.3. For each x G M. d , let p l (dy\x) be the unique invariant measure 
corresponding to the operator C x equipped with periodic boundary conditions in y. 

Under Regime 1, we assume the standard centering condition (see [2]^ for the 
drift term b: 

H x ,y)i J ' 1 (dy\ x ) = 0. 



The variable x is being treated as a parameter here. 

We note that under Conditions 12.11 and 12.31 for each I G {1, . . . , d}, there is a 
unique twice diffcrcntiablc function xi( x iy) that is p~ periodic in every direction 
in y, that solves the following cell problem (for a proof see [2], Theorem 3.3.4): 

( 2 - 3 ) £lxifav) = -h{x,y), \ xi( x ,y)n 1 (dy\x) = 0, l = l,...,d. 



We write x = (Xij ■ • ■ ?Xrf)- With this in hand, it will become useful to define a 
function Xi(x, y), i = 1,2, as follows: 

Definition 2.4. For each one of the Regimes i = 1,2 defined in \2.2\) . let \i : 

R d x y -> M. d be given by 

Ai (x, y) = (I + V v x(x, y)) c(x, y), and 
A 2 (a;,y) = jb(x, y) + c(x, y), 

where x — {xu ■ ■ ■ 1 Xd) is defined by 12. 3\) and I is the identity matrix. 
Moreover, let Xi : M. d — > M. d be given by 

Xi(x) = / Xi{x,y)tf(dy\x), 
Jy 

and let X\{x) be the flow generated by Xi. That is, for each x G M. d , X l s (x) is the 
solution to the ordinary differential equation 

X l t (x)=x+ f Xi(Xi)ds. 
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Note that, due to our assumptions, the ODE for X % is well defined and have 
unique solution in each regime. Moreover, Theorem 2.8 in [5] guarantees weak 
convergence of X e to X 1 in C([0, T]) for any T > 0. Further, it is easy to observe 
that in our case Theorem 2.8 in [5] implies that for all 77 > and i = 1, 2, we have 

(2.4) lim P j sup \X e t - Xt(x )\ > r/1 = 0, T > 0. 



0<t<T 



Our first objective is to understand the limit of the fluctuations process 

Vt= * o e * , ase|0, 

where /3 e is the appropriate normalization rate. Our second objective is to prove a 
limit theorem for an exit problem of X e using the limiting result for the fluctuations 
process. That is, for a smooth C 2 -hypersurface M in M. d , we are interested in 
studying the joint distribution of the hitting time 

t £ = inf {t > : X e (t) g M} , 

and the exit location X e {r e ) g M as e — > under the assumption that r e < 00 with 
probability 1. Precise assumptions on the joint geometry of the vector field A and 
the surface M will be given in Section @] 

We conclude this section with a remark for the degenerate case e/S — >• 0. 

Remark 2.5. In the case e/S — > 0, the results of [5] indicate that the correct pair 
(C x , A(x, y)) is that of Regime 2 with 7 = 0, as long as there is a unique invariant 
measure to the corresponding first order operator. Due to the fact that this operator 
is first order, the existence and uniqueness of an invariant measure is a difficult 
issue and requires additional assumptions on the vector field c(x, y). For this reason 
and for the additional technical difficulties in treating the related Poisson equation 
h3.1\) , we decided not to treat this case in the current paper. See however, Corollary 
\3.2\ for the case 7 = 0, in dimension d = 1 when c(x,y) > 0. 



3. Analysis of fluctuations 

In this section we establish a limit theorem for the correction of X e — X 1 in each 
case. Before stating our results in this direction, we need additional notation. 
Let us consider the auxiliary PDE problem 

(3.1) C % x Ei(x,y) = - (Xi(x,y)-Xi(x)) , / E t (x, y)fi l (dy\x) = 0, 

Jy 

for i = 1,2. Since, by definition, the right hand side of the PDE averages to zero 
with respect to the corresponding invariant measure tf(dy\x), Frcdholm alterna- 
tive implies that the function Si(x,y) is uniquely defined, p— periodic in y, twice 
differentiablc in both derivatives and with bounded derivatives (see Theorem 3.3.4 
in [2]). The function S, will be used to understand the dependence of terms like 

(3.2) 1^ = J* (\i (xt, ?f) - ^(Xt)) ds 
on e and 6. 
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Let us give some preliminary notation. For a function / : R d x y — > M. d , denote 
fi : M. d — > M. d the average with respect to fx 1 : 

/i(af) = / f(x,y)^(dy\x). 

Jy 

Also, define the following functions, 

*i(x,y) = (I + V y x(x,y))$(x,y), 
Ji(x,y) = cV y E 1 (x,y), 

qi(x,y) = {I + W y x)(x,y)cr(x,y)(j T (x,y)(I + W y x) T (x,y), 

and 

9 a {x,y) = (I + V y Z 2 (x,y))*(x,y), 

J 2 (x,y) = (b{I + \7 y Z 2 ) + ^ti-[<j<j T \7 y \7 v Z 2 ]^j (x,y), 

9i{x, y) = (I + V y E 2 )(x, y)a(x, y)<r T {x, y)(I + W y E 2 ) T (x, y). 
Further, for x G R d , let be the linearization of X % along the orbit of x: 

(3.3) i$t(t) = DA*(^)$j.(i), **(0) = z 



where DA Z is the Jacobian matrix of X 1 . We are now ready to state our results. 

leorem 3.1. Let T 

= min ( - 2i sa "l 



Theorem 3.1. Let T > 0, and assume Conditions PH1TOI Set 9{ = ~, 6\ = §-7, 



^ = lim ^— G [0,oo], 



#(4 



R e ,4 = 0, 
\e m ,^€(0,oo]' 

Let 77* oe the Ornstein- Uhlenbeck type of process 

dfjl = DXi{Xi(x ))f] t dt + [^Hiti G (0, 00]) + 1(4 = 0)] Ji(xi(x ))dt + 

+1 (£i ^ 0)) [l (m = ai/2) *j(^*(x ))«ft + 1 (m = 1/2) ^ /2 (X t *(xo))dWi 
% = Col (m = a 2 /2 and 4 ^ 0) . 

Then, for each e > 0, t/iere is a process rf(ii), such that 

Xt=X\+F(li)r(HM) 

holds with probability 1 /or every t > 0, and if {!!■{) — > ff{li), as e — > 0, in distribution 
mC([0,T];M d ). 

In the case 7 = we can give an analogous result for the one dimensional case: 

Corollary 3.2. Let T > and ief the dimension be d = 1. Assume that Condi- 
tion \2.1\ holds, and that c(x, y) > /or every 1 £ 1,1/ G J'. Then, in the case i — 2, 
and 7 = 0, ttie conclusion of Theorem \3.1\ holds by setting in the corresponding 
expressions 7 = 0. 
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The proof of this corollary is omitted since, due to the condition c(x, y) > 0, the 
ODE it = c{x, z t ) has a unique invariant measure for any which then allows 

the proof of Regime 2, presented below, to go through with 7 = 0. 

Before proceeding with the proof of Theorem l3.ll we mention a useful observation 
in the remark below. 

Remark 3.3. Notice that the drift term in the effective equation for fj has an 
extra term Ji(Xl), which is present in the case li ^ oo. This term arises from 
the fluctuations associated with the term 1%' in 13.2(1 . Hence, if li =^ oo, the 
contribution of this term is not negligible in the limit. 

Proof of Theorem \3.1\ Let us first consider Regime 2, i.e., 7 G (0, 00). Let A e (t) = 
XI — X\(xq). so that A e (t) satisfies the equation 



dA € {t) 



X{ 



(3.4) 



+ Veer [Xf, ?±\ dW t , 



,a 1 /2. T; e 



X 



w [Xl,-± -X 2 (X t ) 



dt 



with initial condition A e (0) = e" 2 / 2 £ e . For notational convenience we sometimes 
omit the subscripts i, and xq throughout the proof. 

Note that, since A = A2 is smooth, Taylor's theorem implies that 

A(xi) = X(x 2 ) +D x X(x 2 )(xi -x 2 ) +Q[X](xx,x 2 ), x u x 2 G R d , 

for some function Q[X] such that \x\ — X2\~ 2 Q[X](xi, x 2 ) is locally bounded. Then, 
rewriting (|3.4p . we obtain (with A = A2, A = A2) 



I x 



(3.5) 



b I 

X 

Xi 



xf 

X(X 



ds+ / D x X{X s )A £ (s)ds 



Xi 



ds + e ai/2 / * e [ X 



ds 



41 I <r[Xl-±)dW.+ I Q[X](X s ,Xl)ds. 



To understand the asymptotics of the right hand side of the last display, we need 
to understand the behavior of the integral term 

It = J* X (xi ?f\ - X(Xt) ds. 

For this purpose, apply Ito's formula to "B(x,x/6) = (Bi(x, x/5), . . . , ^(x, x/5)) 
with x = X\. After some algebra, it follows that 



5d~\xir4- 



-11 
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Xf,^l I. 



Xi 



V~e [(V,S + SV x E)a] X t £ , -M dW t + i? e X t e , dt 



x; 



where, for each (x,y) G M. d x y, R e is given by 



R\x,y) 



eW x E + S[c + e ai/2 ^ e ) V X E + etr 



aa T [V x Vy + -V x V x E 



{x,y). 
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Therefore, taking into account the PDE that 3 satisfies, we get that 

it = J' (W 2 vrv y ~ + ey 2 ) (xi, ?f) ds 

- 6 ( a fa f ) " B fa ¥ )) + *f Q ^(x t , f) d w s 

(3.6) +^ R^Xt,^pjd s + 5^ V x Ea(x e s ,^f)dW s 

So from this expression, and the smoothness and boundcdness of the involved func- 
tions we can study use this expression to study the limit of A e . 
Using (|3.6[) in (|3.5j) . and the definition of 1 $>2, it follows that 

A e (t) = e Q2 / 2 f + J D x \(X s )A'(s)ds + J % (x e s , ?f) ds 

+ 6\ J J 2 (xi, ?f) ds + yfeJ o (I + V y S)a(xt, ^f)dW s 

+ jf*fa£)*-<<*-?)- s fa?)) 

(3.7) + S ^J ^^( x h^f)dW s + Q[X\(X s ,Xl)ds. 

We are now going to use this representation to prove that the family {A e / (3 e (£), t S [0, T]} e>0 
is relatively compact in C([0, T]; M. d ). To do so we shall use Proposition 3 of [22], 
which says that the family of processes {?7 e } e>0 is relatively compact in C([0, T]; R d ), 
if there is an eo such that for every h > 0, 

(i) there is a Nh < 00 so that 

P \ sup |?7 t £ | > N h \ < h, e G (0,e ), and 
[te[o,T] ) 

(ii) for every M < 00, 

lim sup P<^ sup \Vt 1 -Vt a \>h, sup ?7 t £ < M > = 0. 

r ^ U ce(0,£ ) [t 1 ,t2&[0,T],\t 1 -t 2 \<r te[o,T] J 

We will prove these two points for the family rf = A e //3 C . 
Let 

®% (t) = e^-™$ X0 (t)C + e^- m $ X0 (t) £ [^(s)}' 1 n (xt, ?f) ds 

(3.8) + ei- m $ X0 (t) f [$ X0 {s)]~ l (I + (x s £ , *f) dW s , 
so that Duhamel's principle implies that 

A e W = e m Q% o {t) + e^ x M f^^MV 1 J2^Xl,^f)ds 

(3-9) +Rim+<S> X0 (t) f [* X0 (s)]- l Q[\](X s ,Xl)ds. 

Jo 



s 
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where the term Rl[&] is defined as 

+ *a (t)£[<Ms)r 1 a e (x^)da 

+ *v^*»o(*) ^ [^o(s)]" 1 V*Sa(xj, ^)dW.. 

Let us show that points (i) and (ii) hold for the family A e //3 e (£), e > 0. First, 
due to boundedness of the involved functions, the definition of R e , and Doob's 
inequality for the martingale term, it is easy to see that 

sup (mvr 1 mm"' 

<t<T \ 



0. 



limE 

c l° lO<t<T 

Likewise, we will show that the family of processes 

A * = WW** ® Jo [^o(s)] _1 Q[A](X,^)ds. 

also converges to uniformly on [0, T] in probability as e — > 0. Once we have these 
two facts, points (i) and (ii) follow for A £ //3 e (£) due to the definition of (3 e {£), 
and the boundedness of all functions. Hence, we are just left to prove that 

sup|A'|^>0, e^O. 

t<T 

Let v <E (1/2, 1) and set 

r» =inf{t: \Xl-X t \ > {f) v }. 
Using the quadratic decay of Q[X], and the fact that 2v > 1, we see that 

lim sup |Aj| = 0, 

with probability 1 . As a consequence, we are left to show that P{T < T e (is)} — > 1. 
To do so, note that (|3.9p implies that if T e {u) < T, 

i = 0ST" sup |A e (t)| 

t<TAr'(u) 

for some random variables C\,Ci < oo, P-a.s. Hence, since the r.h.s of the last 
display converges to 0, it follows that lim e _^o P{ r£ (^) < T} — 0, which implies the 
precompactness of the family {rf — A e /(3 e (£), e > 0}. Clearly, the tightness of the 
family {X e , e > 0} implies tightness of the pair {(r] e , X e ),e > 0}. 

Now that we know that the family {(r] e ,X e ),e > 0} is precompact, we are left 
to identify its limit. Wc shall use p.7[) and the martingale problem formulation. 
An inspection of (|3.7j) shows that the terms in its third and fourth line are of lower 
order compared to the other terms and thus should vanish in the limit. Let us make 
this now rigorous. 

Consider, the process 

q = 4 + 8 (pw)- 1 (b(xi Zf\ -z(x<, ^ 
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and let <j> G C%{R d ). We identify the limit using the martingale problem approach. 
Write for simplicity j3 = f3 e (£). By ltd formula we have 



0(C t £ ) = 0(r 1 e Q2/2 f )+ / D x \(X s )r 
ri 



e 1 



tr 



P Jo 

VV0(C)(/ + V,S)aa T (J + V„S) 
X 



J I XI 



Xi 



xt 



Xi 



ds 



+ 



P Jo 

t r/ e «i/ 2 



V0(C)(/ + V y 3cr) X 



1 



P 



P 



X 



% + -R* )[XI,^ )+-Q[\](X s ,Xt 



dW s 



1 



P 



V0(C)ds+ 



+ 



S 2 el 



tr 



(3.10) 



VV<f>(C)V x E*(V x Zcrr(XI,^- 



Xi 



ds 



It is clear from (|3.7|) that there are different cases to consider, depending on 
the order that the different terms go to zero. For the sake of presentation, we 
shall only study in detail the case m = 1/2 and I ^ 0. In this case the limiting 
process is solution to a stochastic differential equation. The other cases follow 
similarly. Without loss of generality, let us simplify the problem more and assume 
that ai,a 2 > 1. Then, we actually have /3 = P e (i) = \/e and our target is to prove 
that for any < s < t < T 



limE. 



/ [(^A(X)^(r)+r i j(^))v^) 



(3.11) 



+-tr[VV<^)g(X r £ ) 



dr 







This follows directly upon rewriting the left hand side of (|3.1ip using (|3.10p . In 
particular, note that the following holds: 

(i) Due to boundcdncss of 3 



5/3 _1 E sup 

0<t<T 



x >4 



= o. 



(ii) Due to boundedness of the coefficients and the bounds on the derivatives 
of the auxiliary function 3 [2], we have that 



lim E sup 

o<t<rJo 



,«l/2 



P 



xi, 



Xi 



V0(C) 



ds = 0. 



Similarly to the argument used to prove tightness, we also have that 



hmE / (^-^[^(X.X^V^COdr -> 0. 



10 



SERGIO ANGEL ALMADA AND KONSTANTINOS SPILIOPOULOS 



Moreover, it is clear that there is also a constant C < oo such that for the 
remaining Riemann integrals in (13.101) 



E sup 

0<t<T Jo 



ds < C. 



(iii) Due to boundedness of the coefficients and the bounds on the derivatives 
of the auxiliary function S [2], there is a constant C < oo such that, using 
Doob's martingale inequality, for any stochastic integral in (|3.10p 

2 



E sup 

0<t<T 



r 4 



x: 



dW„ 



< C 



(iv) By assumption j3 = s/e and 9/(3 — ?> t 1 as e | 0. 

The validity of p. lip together with tightness of the pair {(77.% Xf), e > 0} and 
uniqueness of the limiting martingale problem imply the claim. 

Let us now concentrate on Regime 1; that is, 7 = 00. Consider the solution 
to the cell problem x = (xi? • • • j Xd), which is periodic in every direction in y and 
satisfies 

£ x Xi{ x ,y) = -b t (x,y), / xi{x,y)V>{dy\x) = 0, I = l,...,d. 



By applying Ito's formula to x{ x i x /$) — {Xi{ x i x /^)i ■ ■ ■ iXd{ x i x /&)) with x = X{, 
we can reduce the problem to the previous case. Indeed, by Ito's formula, the cell 
problem formulation, and Definition 12.41 it follows that 



dx\Xt, 



1 „ e 



tr [aa T V 2 x ] 



;tr [aa T W x W y ] + £ Ql/2 * e • [V x + \v v 



x[xt, 



Xi 



dt 



+ ¥ c 1 xi x(x 



Xi 



dt- 



x[xt, 



A7 



dW t 



Hence, recalling the cell problem (|2.3[) we have 



MX! 



xt 



dt = 



(eb + Sc)-\7 x +c-\7 v + |<5tr [aa T \7 2 x ] 



-etr [aa T V x V y ] + <W 2 * e ■ [V x + -V v 



dt 



+ V~e5 



x[xt,^\dw t -dx[xi^f 



Xi 



Using this in (|3.4p . rearranging terms, and proceeding as in Regime 2, we obtain 
(with A = Ai, A = Ai) 



dA e (t) 



+ V~<I + VyxM Xt, -f )dW t + Fx[ XI 



xt 



x> 



R{(Xt, -f)dW t + Q[\](X tl Xt)dt - 5d x [ XI -± , 



dt 



(i+Vyxmxz, 



x; 
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with A e (0) = e Q2 / 2 £ £ . Here, we have defined 



(eb + 8c) ■ V x \ + etr 



oa L \ -\7 x \7 x + \7 x V y )x 



5e a ^ e (x,y)W x ) 



R\{x,y) = V^SW x x(x,y)cr. 
As in Regime 2, we need to understand the behavior of the term 



-A(X- 



ds. 



For this purpose, apply Ito's formula to S(x,x/#) = (3i(x, x/fi), . . . , S<j(x, x/8)) 
with x = A t c . Taking into account the PDE that 5 satisfies, we get 



x '-f 



ds 



52 '-Jx-M 



\x Q ,— 



where the lower order term R\{t) is 



C2 ft 



Xi 



ds + -e ai/2 I * e V y S( X 



e XI 



ds 



J (bV x Z + ti[aa T V x V y E\) (x e s ,^jds + S 2 ^ J tr [gcj t V x V x E] (x 



+ 5 

5 2 

H 

e 

Hence, we get that 



e / V^Sct X 



Xi 



dW s 



e V y Ea(XI, 



xt 



dW s 



A «( t ) = e a»/a^ + / Dx X{X s )^{s)ds 



X 



ds 



e^' 2 (I + V y x)* e [XI,^)ds + V~e (I + VyXMXt 



X< 



dW s 



ds + I Rl( x l,-±)dW s +Rl(t)+ I Q[X](X s ,Xl)ds 



-II 

6 



x[x 



X 



x ,— 



Using these expressions, tightness follows after using Duhamel principle using the 
same arguments as for Regime 2. The identification and uniqueness of the limiting 
point follows by the martingale problem by considering the process 



Xi 



X! 



X!, 



and we applying Ito formula to a test function cf> £ C 2 (R d ) with stochastic process 
Q. Using the resulting expression, the claim follows by the martingale problem 
formulation as in Regime 2. Thus, the details are omitted. □ 
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4. On Asymptotics for the Exit Time and Exit Location 

Let us state the main result in regards with the correction to the exit. First we 
describe our assumptions on the joint geometry of the vector field A and the surface 
M. We define 

T = inf {t > : X l t (x ) G 3\l} , 
and assume that < T l < oo, for each i = 1,2. Also, for each i = 1,2, we 
denote z l = Xj^(xo) <G M and assume that Xi(z) does not belong to the tangent 
hyperplane T z iM of M at z'. In other words, we assume that the positive orbit of 
Xq under the vector field Ai crosses M non-transversally. 

Let us introduce a local basis around the exit point z 1 in order to express the 
correction. Given z £ M. d , let T Z M be the tangent space of M at z. For i = 1,2, 
define the projections tt 4 : R d -> R and tf M : R d -> T Z M by 

v = 7r*v ■ Xi{z) + irlfV, v e M d . 

That is, 7T l is the (algebraic) projection onto span(Ai(z*)) along T l M and 7r^ f is the 
(geometric) projection onto T z iM along span(Ai(z*)). 

Wc have now all the elements necessary to state the main theorem: 

Theorem 4.1. With the notation of Theorem \4-S\ let £i < oo and assume addi- 
tionally that (~^0l — > 1, for some £ > 0, as e, S 4- 0. Then, in the setting introduced 
above and under these assumptions, for i = 1,2, 

(4-1) ^y(r e - T\X%, - zVKW^V&W)- 

in distribution as e — > 0. 

4.1. Proof of Theorem 14.11 Due to Theorem 14.31 the proof follows similarly to 
the geometric ideas as the proof of Theorem 1 in jTj . Hence, our goal is to show 
that we can apply these ideas to get the proof. However, due to the extra effect of 
averaging, several adjustments are needed. 

Under the assumptions of the theorem, note that f3 e (£) = e@, where 

P = (l(£ = 0) + ml(£^0). 

It follows that in Regime i = 2, f3 = min{C, 1/2, ai/2, 02/2}. 

The following corollary of Theorem 13. II is essential in our proof. After the proof 
of this corollary, we will restate Theorem 13.11 in a more convenient way for the 
propose of this proof. Let us start with the corollary: 

Corollary 4.2. Under the same assumptions as in Theorem \3.1l if, as e — > 0, 

t c -> 0, t e {e/S - 7)e~ m -> 0, and m = a 2 /2, then 

(4.2) mp\e- m (X;-Xi)-?\^*0, e -> 0. 

t<t> 

Proof. First, let us focus on the case i = 2. In this case, observe that, since 
m = a/ 2, from (|3.8[) we have 

where rj is a process such that for any family (s e ) e >o such that s c — > 0, it follows 
that 

sup | r ^(t)|A0, e^0. 
te[o,s'] 
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Hence, from p.9[) it follows that 

+ e- m 6^ xa (t) J* [^(sT 1 (b(I + V„S) + itr [aa T V,V,~]) (x<„ ^ 

+ e- m RI[$] + e- m $ X0 (t) [ [^(s^QiX^XDds, 

Jo 

which implies the result since 9 e t e /e m —> 0, and, as in the proof of Theorem 13.11 
i? e [$] is of order e, and Q[X] of order (/3 € ) 2 as e — s- 0. For Regime 1, the result 
follows analogously from the expression obtained for A c in this case. □ 

We now give a useful restatement of Theorem 13.11 Before doing so, some addi- 
tional notation is needed. For i = 1,2, let 

Ql {t) = l{m = a 2 /2W Xa {t)e 

+ l(m = ai/2)** (t) / K(s)]" 1 *j (XI) ds 
Jo 

(4.3) +l(m=l/2)*i (t) f\$i a {s)]- 1 q] /2 {Xl)dW s , t > 0. 

Jo 

and 

Jo 

We have the following result: 

Theorem 4.3. Assume Regime i = 1,2 and let Conditions \2.lWlK holding. Con- 
sider the stochastic process 

(4.4) nm = et (t)l (4 ^ 0) + flj (t) [irHiU e (0,oo]) + 1(4 = 0)] . 
Then, for each e > ; t/iere is a process rf{&i), such that 

holds a.e. t > wit/i probability 1, and r] e (£i) — > ff(ti), as e — > 0, in distribution in 
C{[0,T};R d ). 

The proof of Theorem 14.31 follows essentially by Theorem 13.11 and Duhamel's 
principle. The details are omitted. 

We are now in good shape to give a direct analog of Lemma 7 in [TJ : 

Lemma 4.4. Let u> £ (/3/2,/3). Then, for each i = 1,2, there are two a.s.- 
continuous stochastic processes T e,± ' 1 such that 

sup |T^|A0, e^O, 

ie[0,e«] 

and a.e. t € [0, e w ] a.s 

(4.5) X^_ t = z l - tliiz*) + eP(ri e T _ t (l) + TJ'"'') 
and 

(4.6) X^ +t = z* + tM*) + (vt(£) + Tt'+' i ) , 
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where fj e (i) is a stochastic process such that for any cu <E (/3/2,{3),£ <E [0, oo], it 
follows that 

(4.7) fj e (£) ff.{£), e ->• 0, m distribution in C([0,5];R d ) for any S > 0. 

Proof. We will omit reference to the regime i when no confusion arises. Let X t (x) ( 
X-t(x)) be the positive (negative) orbit of xq under the flow A. Note, in particular, 
that 

jX i ±t {x)=±\ i {X i ±t {x)l X i (x)=x, 

for every x GM. d . 

Due to Strong Markov property and Lemma |4~51 the process X t e = X^ +T satisfies 
the SDE ([23]) with respect to the Brownian Motion W(t) = W(t + T) - W{T) and 
initial condition Xq = z + e^r]^,(£). Hence, we can apply Lemma |4~51 to this shifted 
equation, to obtain 

where /? = £ if i = 2, rh = min{/3, 1/2, ai/2} > £, and $ = rh in other case. The 
process rf is such that rf{t) ->~t; 4 (£), as e — > 0, in distribution in C([0, S}; R d ) for 
any S > 0, where rf is defined by 

(4.8) fji(£) = ei(t)i (£^o) + [i-Hii € (0, oo)) + 1(£ = 0)] H l z (t), 

where 

ei(t) - K(t)vW) 

+ l(m = ai/2)*t(f) / [$t(s)] _1 *(Xi(2))d S 
Jo 

+ l(m = 1/2)$* (i) /' [^(s)] _1 9- /2 (X:( Z ))dW s , 
Jo 

and £ = lim e _»o e m (e/S — 7) . Note that £ = if and only if m > f. In this case, 
/3 = C < m < min{l/2, ai/2}. Hence, in any case, rh — f3, and, from the definition 

of ft j8 = /?■ 

Using these results, it follows that 

X$. +t = z l + t\(z l ) + e?{r,t{£) + T^) , 

where T t e ' + ' ?; = e - /3 (X*(z 4 ) - (z 1 + tX(z 1 ))) . Hence to show that gU) follows, we 
have to prove the convergence of T e ' + ' 1 towards zero and (|4. 7[) . Let us start with 
the former, to show that T e,+ ' 1 converges towards 0, use Lemma 6 in p] to find two 
positive constants C\ and C2 such that for any t > 0, and x <E M. d , 

(4.9) sup \X* ±t (x) - (z' l ±tX(z))\ < C ie C2t (t\x- z l \+t 2 ). 

s<t 

This implies the convergence to of Y e ' +,? . We are left to prove (14.71) . In order to 
so, recall that m = /3, and Corollary 14.21 implies (|4.7p . 
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To get P~5]l , note that, for t <G [0,T], Lemma B~31 Equation (g^J), and the fact 
that X^,_ t (xo) = X z _ t (z l ) imply that 

X^_ t = X^ T _ t (x ) + e^ T _ t (e) 

= X i _ t (z i )+^V e T-t(t) 

where T^' - '' = e"^(Xi t (^) - (V - tX^z 1 ))) converges to due to (@~9]). □ 

From this lemma, the proof follows using similar arguments as in pQ . We start by 
giving an alternative characterization of the time t £ that will be used to explicitly 
write out the correction X^. e — z % . 

Parameterize the hypersurface M as a graph of a C 2 -function F % over T Z M, i.e., 
y M> z l + y + F l (y) ■ \i(z) gives a C 2 -paramctrization of a neighborhood of z % in M 
by a neighborhood of in T Z M. Moreover, DF(0) = so that \F(y)\ = 0{\y\ 2 ), 
y —> 0. With this definition, it is clear that, for w 6 M d with w — z small enough, 
w £ M if and only if ni,(w — z) = F(ttm(w — z)). Moreover, Lemma l4.4l and a direct 
application of Lemma 8 from p] gives that for every w € (/3/2,/3), 

liniP({r^f e }n{|r e -r|<e"})=l, 

where f e = inf{i > : tt 1 (X* - z l ) = F (tt 1 m (X{ - z 1 ))}}. Hence, from now on, 
we fix ui <E (/3/2,/3) and condition on the intersection of the events Slf = {f e = r e } 
and 0| = {|r £ -T 4 | < e w }. 

Let f e — T l ~ t c , so that r e — T l ~ f e , and (since we are conditioning in fi| ) 
|f e | < e". In case f 6 > 0, the conditioning on Of , and (|4.5[) imply that 

~r + (p e Ti _ f . (£) + T'T'*) = , 

where D e '+ = —f^A^z) + e^ 71 !/ (^-f « W + T f=~' 1 )- Hence, since 2uj > /3, it 
follows that 

(4.10) l({r l >TjnO e 1 n^)e~ /3 ((T e -T l )-7r 4 77^_ fe (£)) Ao, e^O. 

In the case f e < 0, the reasoning is completely analogous. Indeed, using (|4.6p . we 
can get that 

(4.11) l({T 4 <TjnO e 1 n^)e^' 3 ((T e -T 4 )-7r l fyl fe (^)) A 0, e -> 0. 
By adding up (|4TTU|) .and P~TTj) . and using (|4.7[) in (|4. 10[) it follows that 

e -/3(V - T l ) nhf T (£), in distribution, as e 0. 

Once we have the time component, the spatial component and the joint conver- 
gence follows as in [T]. 

5. First Order Langevin Equation in a Rough Potential 

In this section we apply Theorems 13.11 and 15.31 to a small noise diffusion process 
in a rough environment in Regime 1, i.e. we assume that <5 goes to zero faster than 
e does. We change the notation to include the dependence on 5. 
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To be precise, we consider the first order Langevin equation in a rough potential, 
defined as 



(5.1) dX e / 




W ( X e / 



dt+V2eDdW t , X^ 



Xn. 



where Xq — > x as e 1 0. Note that the rough potential V e ' 5 (x, f ) = eQ(x/S) + 
V(x) is composed by a large scale smooth part , V(x), and by a smaller scale fast 
oscillating part, eQ(x/S). We assume that V is a C 3 (IR d ) function while Q is C 2 (M d ) 
with period p. 

By (|2.4p we know that X e - S converges in probability, uniformly in t <G [0,T], as 
e/5 j" oo, to X where X is the solution to the ODE 

X t — A (X t ) , X = XQ, 

driven by the vector field A defined in Definition 12.41 In this case is easy to see 
that the invariant measure n(dy) of C = —VQd x + Dd x is given by the Gibbs 
distribution fj.(dy) = K~ 1 e~~^ L dy, where 

K = I e D dy. 
Jy 

In dimension d = 1 , after some algebra, we get that 

A(X) = — CrV'(x), 

KK 

where K = f^e^v^dy. 

Let us first see how the central limit type of Theorcm l3 . 1 I translates in this special 
case of interest. In many problems of interest, one is interested in understanding 
the behavior of the process starting within the neighborhood of a stable point of 
V^(x), assume that such a point is x = xo. In Figure [I] we see a simple example of 
such a potential function. To account for this fact we assume that Xq = x + e° 2 / 2 £ e , 
where the random variable £ e — > £° in distribution as e \. 0. 

The following proposition states the central limit theorem in this particular case. 
We point out the presence of the additional drift term Ji(X t ). 

Proposition 5.1. Consider the solution to the SDE \5.1\) in t <G [0, T] where 
e/5 j" oo. Under the setup of Theorem \3.1\ we have that the process = ^e^* 

converges in distribution in the space of continuous functions inC ([0,T];M d J to the 
process fj t which is as follows. 

(i) If I = 0, then fj satisfies the ODE 

dfjt = [DX(X t )fj t + Ji(X t )] dt, fj = 

(ii) If £ € (0, oo], then fj is solution to the Ornstein-Uhlenbeck process 
dfj t = [DX(X t )fj t + i^MXt)] dt + l{m = l/2}q{ /2 dW t , fjo = £°l{m = a 2 /2}. 

In dimension d = 1, we have that A(x) = ~~^~gV\ x ) an d Qi = — P ^P and 
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-2-10 1 2 

X 



Figure 1. Q{x/8) = cos{x/5) + sm(x/S), V(x) = \x 2 with e = 
0.1,<5 = 0.01. 

The proof is a straightforward application of Theorem 13.11 and thus omitted. 
Notice that unless Q = or I = oo, the term J\{x) ^= has non zero contribution 
in the limiting fluctuation process. 

Next, we study the related conditional exit law, the result is in Theorem 15.31 
From now on we assume that the initial point is Xq = xq and restrict the analysis 
to dimension d = 1 and to taking first S \. with e fixed and then taking e J, 0. 
Essentially, this corresponds to the case i = oo. 

Let us assume for concreteness that V{x) is strictly convex, has unique minimum 
at x = z such that V(z ) = V'(zo) = 0, V(x) > for x ^ z and V'(x) ^ for 
x =/= zq. Without loss of generality we assume that zq = 0. Consider an interval 
/ = [it_,.t + ] containing xq and assume that < X—. In Figure [lj we see a simple 
example of such a potential function. 

Let us define the exit time 

t c - s = inf {t > : X e ' S {t) i {x-,x + )} , 

and consider the event B f - S = ^X*'^ s = .t_|. From the assumptions of V, it follows 

that lim e ^o P {B e ' S } = 0. Large deviations for such and more general processes of 
similar structure have been studied in [5]. Our goal in this section is to study 
the behavior of the exit problem from / for this process conditioned on the rare 
event B e ' S when S goes to zero much faster than e and provide a limit theorem 
using Theorem 14.11 With some abuse of notation, we shall denote this process as 
X £ < s \ B e,5. In Remark El we discuss the case when Q is a stationary and ergodic 
random field. 

We investigate how the fast oscillations of the small perturbation function Q(y) 
affect the conditional exit law. We first derive the process to which X e ' S \ Be ,s con- 
verges to by first taking 5^0 with e fixed and then taking e j 0. We do this 
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using the Feller characterization of one-dimensional diffusion processes given in [5] 
and the related weak convergence results of [11]. These results are recalled for the 
convenience of the reader in Appendix [X] 

In order to formulate our results we need to introduce some more notation. Let 
us define the 1-dimensional torus that the fast motion takes place as T 1 . We shall 
use the notation, 

(5.2) (g) = - [ g(z)dz 

P Jti 

for the mean value of a periodic function g with period p. 

Theorem 5.2. Let S > be given. Given the assumptions made on the functions 
Q(y) and V(x) above, we have that conditioned on B e ' S , the process X e ' S converges 
weakly, as S 1 with e fixed, in the space of continuous function C ([0,5];R) to a 
process which at least up to the time that it exits the interval I = x+], satisfies 



(5.3) dXl 



V 



'(*) 



dt 



'2D 



-3. 
e o 



e i> 



--dW t , X Q 



x 



where the function ty e (x) = o(e) as e X uniformly in x. 

The proof of this theorem is deferred to the end of this section. With this result 
at hand, Theorem I4.1[ implies the following result for the limiting distribution of 
the conditional exit time r e ' S . 



Theorem 5.3. Let 



T(x ) 



dy 



< oo. 



V'(y) 

Given the assumptions made on the functions Q(y) and V{x) above, we have that 
conditioned on B t S , the distribution of -^=(r e ' S — T(xq)) converges weakly, as first 
5^0 and then e J, to a Gaussian random variablee with mean zero and variance 
given by 



2D 



Q 

e ° 



1 



-dz. 



Proof. It follows along the lines of Theorem 2 in pQ and thus omitted. 



□ 



The results hold in the case of a random environment. In particular we have the 
following remark. 



Remark 5.4. Even though we have stated Theorems \ 5.2\ and \5.3\ only for a pe- 
riodic function Q(y), the proof of Theorem below immediately shows that the 
statements are true also when Q(y) is a stationary, ergodic random field defined 
on some probability space (ty,Q,i/). For every ui € Q(y,uj) is C 2 (IR) in y with 
bounded and Lipschitz continuous derivatives up to order 2. In particular, in this 
case we have 



= E v 



Q(y) 



= E v 



where E u is expectation under the random environment. In the case of Theorem 
\5.S\. it is easy to see that the convergence is weak in C ([0,T];]R) ; in probability with 
respect to v. 
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From Theorems 15.21 and 15.31 we can get some interesting conclusions on the 
effect of the small but fast oscillations eQ(x/5) on the underlying potential V(x). 
We have the following remark 

Remark 5.5. Theorem \5.3\ gives a second order approximation for r e ' S conditioned 
on B e ' S when e, S are small. We notice that compared to the process without any 
fast oscillations (i.e. take Q(y) = 0), the standardized limiting conditional exit law 



has variance multiplied by the constant 
is easy to see that 



Q 

e d 



By Holder's inequality, it 



> 1 



Therefore, the limiting conditional variance has been enhanced by the factor 
due to the fast oscillations. 



_Q 

e d 



We conclude this section with the proof of Theorem 15.2 



Proof of Theorem ] 5. 2i By Lemma 3 in [T] we know that conditioned on B e ' S , X e ' 5 
behaves, for each e, S > 0, as a diffusion process with infinitesimal generator L £ ' S 
given by 



where 



and 



(*!)--!<>' (j)- **<«>+ *° if 



K) 



d x + eDdl 



h e > b {x) 



Jo h^(y)dy 



h e ' s (x) =exp 



1 

~DJ Q 



--Q' 



V'(y) 



dy 



Then, as it is ease to see, the operator L e,s can be equivalently written in the 
D v e,sD u t,s characterization of Feller [5] (see Appendix 1X1 for some related results 
from the literature). In this case, the corresponding u e ' S (x) and v e ' S (x) functions 
are defined as 



u e < s (x) 
v e ' s (x) 



[ cxpi / b e ' S (z, -)dz \ dy, and 

Jo I £ D Jo 6 J 



o 



cxp 



1 

7b 



¥> s {z, ^)dz \ dy 



By [H], we know that if u e (x),v e (x) are the limits of u e,s (x), v e ' S (x) as S I 0, 
then the process corresponding to the operator D v c,sD u c,s will converge weakly in 
C ([0, S}; R) to the process corresponding to the operator D v zD u i. Our task now is 
to investigate these limits. 



20 SERGIO ANGEL ALMADA AND KONSTANTINOS SPILIOPOULOS 

Let us first investigate u € ' S (x). To simplify notation, denote = — — jp ; 

^(x) = — j and £(?/) = <fr(p)dp. With this notation at hand, we have 

=/"«p{-i/"^ I J)cfc}* 

CX P(-C(?)-- /%(*)<**) rap (-2 r^r^L rf 4^ 



o J Q z h^(w)dw 



V 



-c(f) - 7 / 



By the mean value theorem we know that for a periodic function g G L a (y), for 
a > 1, we have that .9 (f ) — (g) in L[ a oc (R) as <5 i 0. The convergence is in the weak 
sense in the spaces of functions for any arbitrary bounded interval in R. Using this 
on (|5.4|) . we have that as S i 0, u e '*(a;) — > u e (x), where 

o f* (e-C)e-iK- i '^ d Pdw dZ f ^ 



u e (x) = 


J (e c )cxp<i 


' 1 


10 


dz j 


(5.5) 


J (e _c )exp< 


'_1 


10 





Due to our assumptions, we have that J Q y \f (2) <iz < for y > 0. Thus, it is easy 
to see that as e 1 

In (e"* J? - l) = -1 /% (z) dz + ^ e (y) 

where ^ rC (y) = o(e) uniformly in y. Hence, this and (|5.5[) implies that as e J, we 
have 









. Jo 










( exp< 














[-K- 







On the other hand, in exactly the same way we get that v e ' S (x) — > v e (x), 5 J, 0, 
where 



*(z)dz + * e (y) ^ 



with the same \£ e (j/) = o(e) as e | 0. 

The limiting it e (x) and w e (a;) correspond to a process characterized by the gen- 
erator (see Remark I A. Il l 

L e = lb(x) + y f (x))d T + ead'i 



(b(x) + * e (x)) d x 
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where 

2D 

I _Q\ I Q 

e d e" 



b{x) = 



-Q\ / Q 



which concludes the proof of the theorem. □ 



Appendix A. Feller's Characterization of One Dimensional 

Diffusions 

Consider a stochastic process in one-dimensional characterized by its generator 
(A.1) ^) = ^(x)g + 6(x)| 

with smooth enough coefficients a(x) > and b(x). For the convenience of the 
reader, we briefly recall the Feller characterization of all one-dimensional Markov 
processes, that are continuous with probability one (for more details see [8]; also 
|21j). All one-dimensional strong Markov processes that are continuous with prob- 
ability one, can be characterized (under some minimal regularity conditions) by a 
generalized second order differential operator D v D u f with respect to two increas- 
ing functions u{x) and v(x) and its domain of definition. In particular, u(x) is 
continuous and v(x) is right continuous. In addition, D u , D v arc differentiation 
operators with respect to u(x) and v(x) respectively, which are defined as follows: 
D u f(x) exists if D+f(x) = D~f(x), where the left derivative of / with respect 
to u is defined as follows: 

fix - h) - fix) 

£), fix) = lim '— ; j— r provided the limit exists. 

HO u(x - h) - u(x) 

The right derivative f(x) is defined similarly. If v is discontinuous at y then 

Jyyj Ho v(y + h) - v(y - h) 

Remark A.l. For example, it is easy to see that the operator L in II A . 1\) can be 

written as a D V D U operator with u and v as follows: 

_ fy 2H±l dz I ~ 2 rv B 
2 J o a ( z ) nil ana in r = / p J o 



(A.2) uix) = / e J ° dy and vix) = / ——re J ° ^ dy. 

Jo Jo a (y) 

The representation of u(x) and v(x) in HA . ty) is unique up to multiplicative and 
additive constants. One can multiply one of these functions by some constant and 
divide the other function by the same constant or add a constant to either of them. 

Another useful result in this direction is that of , where it is proven that if we 
have a sequence of operators {D v <,D ut ,e > 0} uniquely characterizing a sequence 
of Markov processes {X e , e > 0} such that u e (x) — > u(x) and v e (x) — > v(x) as e 4- 0, 
such that D V D U corresponds to a strongly continuous homogeneous Markov process 
X, then X e X as e I in distribution in C([0, T];R) for every T > 0. 



22 



SERGIO ANGEL ALMADA AND KONSTANTINOS SPILIOPOULOS 



References 

[I] Sergio Angel Almada, Yuri Bakhtin, Scaling limit for the diffusion problem in the Levinson 
case, Stochastic Processes and their Applications, (2011). 

[2] A. Bcnsoussan, J.L. Lions, G. Papanicolaou, Asymptotic Analysis for Periodic Structures, Vol 
5, Studies in Mathematics and its Applications, North-Holland Publishing Co., Amsterdam, 
1978. 

[3] D. Baier and M. I. Freidlin, Theorems on large deviations and stability for random perturba- 
tions, Dokl. Akad. Nauk SSSR 235 (1977), 253-256. = Soviet Math. Dokl. 18 (1977), 905-909. 

[4] J. D. Bryngelson, J. N. Onuchic, N. D. Socci, P. G. Wolynes, Funnels, pathways and the energy 
landscape of protein folding: A synthesis, Proteins, Vol. 21, Issue 3, (1995), pp. 167-195. 

[5] P. Dupuis and K. Spiliopoulos, Large deviations for multiscalc problems via weak convergence 
methods, Stochastic Processes and their Applications, 122, (2012), pp. 1947-1987. 

[6] P. Dupuis, K. Spiliopoulos, H. Wang. Importance sampling for multiscalc diffusions, SIAM 
Journal on Multiscale Modeling and Simulation, Vol. 12, No. 1, (2012), pp. 1-27. 

[7] P. Dupuis, K. Spiliopoulos, H. Wang. Rare Event Simulation in Rough Energy Landscapes. 
2011 Winter Simulation Conference, appeared. 

[8] W. Feller, Generalized second-order differential operators and their lateral conditions, Illinois 
Journal of Math. 1, (1957), pp. 459-504. 

[9] M.I. Freidlin, The Averaging Principle and Theorems on Large Deviations, Russian Mathe- 
matical Surveys 33:5, 1978, 117-176. 

[10] M.I. Freidlin, R. Sowers, A comparison of homogenization and large deviations, with appli- 
cations to wavefront propagation , Stochastic Process and Their Applications, Vol. 82, Issue 
1, (1999), pp. 23-52. 

[II] M.I. Freidlin, A.D. Wentzell, Necessary and sufficient conditions for weak convergence of 
one-dimensional Markov process, The Dynkin Festschrift: Markov Processes and their Appli- 
cations, Birkhauscr, (1994), pp. 95-109. 

[12] M.I. Freidlin, A.D. Wentzell, Random Perturbations of Dynamical Systems, 2nd Edition, 
1998, Springer. 

[13] A. Guillin, Averaging principle of SDE with small diffusion: moderate deviations, Annals of 

Probability, Vol. 31, No. 1, (2003), pp. 413443. 
[14] R.Z. Hasminskii, Stochastic Stability of Differential Equations, Sijthoff and Noorhoff, 1980. 
[15] C. Hyeon, D. Thirumalai, Can energy landscapes roughness of proteins and RNA be measured 

by using mechanical unfloding experiments?, Proc. Natl. Acad. Sci. USA, Vol. 100, No. 18, 

(2003), pp. 10249-10253. 
[16] W. Janke, Rugged Free-Energy Landscapes, Lecture Notes in Physics, Volume 736/2008, 

Springer, (2008). 

[17] Klebaner, F. C. and Liptser, R. Moderate Deviations for Randomly Perturbed Dynamical 
Systems, Stochastic Processes and their Applications, v. 80, 1999, pp. 157-176. 

[18] S. Lifson, J.L. Jackson, On the self-diffusion of ions in a polyelectrolyte solution, Journal of 
Chemical Physics, Vol. 36, (1962), pp. 2410-2414. 

[19] Liptser, R. Sh. and Shiryayev, A. N., Theory of martingales, Mathematics and its Applications 
(Soviet Series), 49, 1989. 

[20] Liptser, Robert and Stoyanov, Jordan, Stochastic version of the averaging principle for dif- 
fusion type processes, Stochastics Stochastics Rep., Stochastics and Stochastics Reports, 32, 
1990, 3-4,145-163. 

[21] P. Mandl, Analytical treatment of one- dimensional Markov processes, Springer: Prague, 
Acadcmia, (1968). 

[22] E. Pardoux, A.Yu. Veretennikov, On Poisson equation and diffusion approximation 2, Annals 

of Probability, Vol. 31, No. 3, (2003), pp. 1166-1192. 
[23] G.A. Pavliotis, A.M. Stuart, Multiscale methods: Averaging and Homogenization, Springer, 

2007. 

[24] J.G. Saven, J. Wang, P.G. Wolynes, Kinetics of protein folding: The dynamics of globally 
connected rough energy landscapes with biases, Journal of Chemical Physics, Vol. 101, No. 
12, (1994), pp. 11037-11043. 

[25] R. Zwanzig, Diffusion in a rough potential, Proc. Natl. Acad. Sci. USA, Vol. 85, (1988), pp. 
2029-2030. 



SCALING LIMITS AND EXIT LAW FOR MULTISCALE DIFFUSIONS 



23 



Department of Statistics and Operations Research, University of North Carolina, 
304 Hanes Hall CB #3260, Chapel Hill, NC 27599 
E-mail address: salmada3@unc.edu 

Department of Mathematics & Statistics, Boston University, 111 Cummington Street, 
Boston MA 02215 

E-mail address: kspiliopSmath.bu.edu 



